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Abstract 

This paper can be thought of as a remark of [5], where the authors 
studied the eigenvalue distribution fix N of random block Toeplitz band 
matrices with given block order m. In this note we will give explicit 
density functions of lim fJ.x N when the bandwidth grows slowly. In 

N— too 

fact, these densities are exactly the normalized one-point correlation func- 
tions of m x m Gaussian unitary ensemble (GUE for short). The series 
{ lim ^x N \m G N} can be seen as a transition from the standard normal 

iV— >oo 

distribution to semicircle distribution. We also show a similar relationship 
between GOE and block Toeplitz band matrices with symmetric blocks. 

Keywords: block Toeplitz matrix, GUE, GOE, limit spectral distri- 
bution 



1 Introduction 

A block Toeplitz matrix is a block matrix which can be written as 



where {A_^-i)> Ao, Am-x\ are m x m matrices and A s = (ay (s))™- =1 . 
If the a,j(s)'s are random variables, then we call TV a random block Toeplitz 
matrix. We suppose the ay(s)'s are real random variables and A s — (A^ S ) T . 
Besides, we list some assumptions as follows. 
Independence of the elements: 

(1) ai 1 j 1 (si) and ai 2 j 2 (s2) are independent if |si| ^ | | , 
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(2) If s 7^ and («i , Ji ) ^ (}2, 32) then ai 1 j 1 (s) and aj 2j - 2 (s) are independent, 

(3) If (ii.ii) 7^ («2,J2) and (ii.ji) ^ (.72, £2) then a ili:t (0) and a i2n (0) are 
independent. 

Uniform boundness condition: 

(4) E[ay(s)] = 0, E[|a y -(s)| 2 ] = 1, -(JV - 1) < s < N - 1, 1 < i,j < m 
and 

sup {E[|a i:) (s)| fe ]|l < i,j < m} = C k . m < +00 

-(N-l)<s<N-l 

and 

sup C k , m — C k < +00. 
Slowly growing bandwidth condition: 

(5) The block Toeplitz matrix is a band block matrix with bandwidth bpf, that 
is, A s = for \s\ > &jv- Moreover, 6jv satisfies: lim 6 at = 00 and 6tv = o{N). 

N—±oo 

Gaussian unitary ensemble (GUE for short) (W m , dfx) is the space T-L m of all 
Hermitian (m x m)-matrices with a certain Gaussian measure d/i (see (TJ [6]). 
We will regularly use the notations (-)gue and (-)tbm, respectively denoting 
the expectations under GUE and random Toeplitz band matrices. If H = (hij) 
is a matrix from GUE, then it is easy to see (/ijj/ij^GUE — 1 while all other 
second moments are equal to zero: (hijhki}GXjE, = 0, whenever ^ (l,k). 

If A — (a,ij (oj))fj = i is a real symmetric or complex Hermitian random matrix 
and its entries are random variables on a probability space il with a probability 
measure P, then the eigenvalue distribution of A is 




where Aj(w)'s are the N real eigenvalues of A. We have 

Theorem 1.1. Let Tjv be an mN x mN random block Toeplitz matrix as 
above. Set Xn = ' ^ len M-Yn converges weakly to f m (x)dx as N — > 00. 

Moreover, if the bandwidth satisfies X}jv=1 ^J? < °°i then the convergence is 

m — l 

almost sure. Here f m (x) — -j= tp 2 {^/mx) where ipj is the jth nomarl- 

Vm j=o 

-4- 

ized oscillator wave-function: ipj( x ) = ~ — 7 J arl ^ ^3 * s ^ e J^ 1 Hermite 

V \/27rj"! 

polynomial:Hj(x) = (— lpe^ ^fr e_ ^~ • 

Remark. f m {x) is the one-point correlation function of GUE up to the scaling [BJ. 
As /1 is the density function of the standard normal distribution and when m — > 
00 fm converges in law to the semicircle distribution (see [5J), { lim iix N \m € 

N—too 

N} can be seen as a transition from N(0, 1) to semicircle distribution. 

In [4], Kologlu, Kopp and Miller got a very similar result. They proved that 
the limiting eigenvalue density of symmetric block circulant Toeplitz ensemble 
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is the same as the eigenvalue density of GUE (Theorem 1.4(1) of [3]). But 
there is an essential difference between their model and our model: In [4] the 
block Toeplitz matrix has a circulant structure while in this paper the block 
Toeplitz matrix has a band structure. In the viewpoint of combinatorics, both 
of the two structures are strong. So they are two different models induced 
from the block Toeplitz matrix model and interestingly have the same limiting 
eigenvalue distribution. There is also another difference: In [4] the entries of the 
block Toeplitz matrix are i.i.d random variables but in our model the entries 
only have to satisfy the independent condition. 

Theorem 1.2. Let v\, . . . , v r be nonnegetive integers, and v = v\ + • • • + v r . 
Set Y N = -5^, then 

r r 

lirn^ _(J](trFA)^>TBM = (Q^Dgue. (1.1) 

°° i=l i=l 



2 Proof of Main Results 

Definition 2.1. Let [n] = {1, 2, ra}, Vn G N. 

r 

(1) We call 7r = {{ai, b±}, {a,k, &fc}} a pair partition of [2k] if (J {a,, bj} = [2k] 

and {ai, b{\ n{ a i? bj} = if i ^ j . For convenience, we assume that a± < • • • < 
at and a, < bi for each pair, under which such a pair partition is called a Wick 
coupling. For 7T, we define ir(a,j) = bj and ir(bj) — cij (1 < j < k). We denote 
by Vi{2k) the set of pair partitions of [2k]. 

(2) Suppose 7r G ViQk). Then 7r can be seen as a permutation: (ai, b\) ■ ■ ■ (a^, 6^) 
Consider the canonical cycle 70 = (1, 2, 2fc— 1, 2fc). Let g(7r) denote the num- 
ber of orbits of the permutation 70 o 7r. 

The following lemma is a well-known consequence of Wick's formula on the 
moments of GUE. One can get it with the method of moment generating func- 
tion (see Section 3.3.1 of [1]). 

Lemma 2.2. Suppose H — (6y)7j'=i is an m x m random Hermitian matrix 
from GUE. Set Y = -j=- Then the odd moments of fly ore all and the even 
moments of [ly are ni2k(^Y) — X) m 9 ^'^ 1 . 

TreV 2 (2k) 

Proof of Theorem ] 1. 11 From Theorem 4.3 of [5] and Lemma 12.21 we know that 
/iy and 7™ have the same moments. It follows from Carleman's Theorem (see 
[3]) that 7™ should be fiy- So the density function of 7™ is the one-point 
correlation function of GUE, which is a well-known function [HE]: 

^ m — 1 

f m (x) = —= V ipj(Vmx). 
From [2] we know if X)iv=i < 00 > then the convergence is almost sure. □ 
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Proof of Theorem ] By Lemma 2.2 of [5], the main contribution in the ex- 
pansion of ■^(n[=i( tr ^ / Af)' yi )TBM comes from the pair partitions. As in the 
proof of Theorem 4.3 in [5], we can show that 







if i s odd 



JV-s-oo N 



m f '"' if ^ ivi is even 

r i=l 



Here the definition of F(tt) as follows: Suppose that ivi = 2v is even and 
7r = {{ai, b%}, . . . , {a„, bp}} € Vsftp), then i^vr) denotes the number of "free in- 
dices" of the system < a * ^ bl ' > (1 < i < £>), where / is defined as below: 

r 

Set i^o = 0. For 1 < x < i^j, if 3 < s < r — 1 and 1 < a < i^ s +i such 



that a; = ivi + (s + l)a, then f(x) = ^ ivi + (s + l)(a — 1) + 1; otherwise 

i=0 i=0 

f(x)=x + l. 

By Wick's formula, we can compute the integral (01=1 (trif'^^GUE an d 
complete the proof. 

□ 



3 Similar Results for GOE and Block Toeplitz 
Matrix with Symmetric Blocks 

First, we remark that we can get the same results associated with GUE if one 
of the following conditions are imposed: 

1) each block of Tn is a complex matrix and A- s = (A S ) T ; 

2) each block of Tn is a Hermitian matrix, that is, A- s = A s = (^4 S ) T - 
However, the situation becomes different when each block of Tn = iAi—i)i t j=i 

is a symmetric matrix, that is, A_ s = A s = (^4 S ) T - Besides, we modify the sec- 

{1 VtiJzj 
, Vs. 
2 if i = j 

We still use the notations (-)goe and (-)tbm, respectively denoting the ex- 
pectations under GOE and Toeplitz band matrices with symmetric blocks. 

Theorem 3.1. Let Tn be an mN x mN random block Toeplitz matrix as men- 
tioned above. Set Xn = ,„ Tjv , , then 

lim _L(trX£} TB M = -(tr(ff/V^) fc ) GOE . (3.2) 
N^roo ml\ m 
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Moreover, fix N converges weakly to g m {x)dx as N — > oo. And if the bandwidth 
satisfies Ejv=i ^J? < °°i fieri the convergence is almost sure. Here 

1 m ~ poo 

g m (x) = —=y^tf(y/mx)+( — ) 1/2 ifj m ^ 1 (y/mx) / s(x-t)i/j m (y/mt)dt+a m (x), 
v 3=0 J oo 

(3-3) 

e{x) — ^sign(x) and ipj is the jth nomarlized oscillator wave-function as in 
Theorem \l.l\ while 



C^m (^) 



±ip 2s (\/mx) -j- ip 2s (Vrnt)dt if m = 2s + 1 
if m = 2s 



Proof of Theorem\3l\ For iJ = let y = ^L. From Wick's formula 

we know ^(tr(iJ/ x /ro) fe ) G OE = when k is odd and ~(tr(H / 'y/m) 2k ) GO E is 



m * 



1 X! ^ {(ht a ,t a , +x ht^t h , +1 )GOV ■ ■ ■ {ht ai .t ai . +1 ht bl ,t w+1 ) gov) 



*i,...,t 2fc =i 



where the second sum is taken over all n = {{a l5 bi}, {a^, G V2(2k) and 
*2fc+i := *i- is symmetric, so {h tai t* 1+1 h tbl t H+1 ) goe ■ • • (^t^+i^t^+^GOE 7^ 

implies < a * 6, ~'" 1 or < bi for 1 < i < fc, moreover, each term 

/, , \ I 2 if t 0i = tbi = ioi+l = ibi+i 

11 II otherwise 

For a given 7r = {{oi, h}, {a fc , fofc}} G V 2 (2k), set 

A(tt) = {(ti,..., t 2fc ) G [m] 2k j^=^+ 1 or{ <a ' = ^ for 1 < i < k\. 



For t = (ti,...,t 2 k) G A(tt), let r(7r,t) = fl{i G [fc]|t„ 4 = i 6j = £ 0i+ i = 
then i<tr(ff/V^) 2fe } GO E = m-^ 1 £ £ 2^'*). 

TTG-P 2 (2fc) teA(Tr) 

Similarly as in the proof of Theorem 3.2 and Theorem 4.3 of [5], we know 
when k is odd ^(trX^) T BM = o(l) and ^(trX^ fc ) TB M is 

2-2. 2. (2mb N )*.mN 11 W l + 2^0 g J, 

iv 2 r ( ,r - t ) b ™ 2fe ' 

= E E E (2mbjv)fc . mjV E nw*+E e ^Kw)+ o « 

OT"(7r,t) 

-EE W^ 00 )- 

TreV 2 (2k) t6A(7r) 
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Thus lim (trX^) T BM = ±(tr(H/^m) k } GO E- From [2] we know the con- 

N—too 

vergence is almost sure if X)jv=i < 00 • Finally, g m { x ) is the one-point 
correlation function of GOE with order m, thus we complete the proof. 

□ 

Remark. The density (|3.3[) follows from the 1-point correlation function of GOE 
(see (7.2.32) of [B]). This family densities can also be seen as a transition from 
the normal distribution N(0,2) to the semicircle distribution with variance 1. 
In the situation of GOE, we also have parallel results to Theorem 11.21 
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